PROCEEDINGS 

OF SCIENCE 



Center Vortex Model for SU(3) Yang-Mills Theory 



Markus Quandt , Hugo Reinhardt 

Institut fur Theoretische Physik 
Universitat Tubingen 
E-mail: |quandt@tphys . physik .uni-tuebingen.de 



Michael Engelhardt 

Physics Department 

New Mexico State University, Las Cruces 



E-Mail: engel@nmsu.edu 



The center vortex model for the infrared sector of SU (3) Yang-Mills theory is reviewed. After 
discussing the physical foundations underlying the model, some technical aspects of its realisation 
are discussed. The confining properties of the model are presented in some detail and compared 
to known results from full lattice Yang-Mills theory. Particular emphasis is put on the new phe- 
nomenon of vortex branching, which is instrumental in establishing first order behaviour of the 
SU (3) phase transition. Finally, the vortex free energy is verified to furnish an order parameter 
for the deconfinement phase transition. It is shown to exhibit a weak discontinuity at the critical 
temperature, in agreement with predictions from lattice gauge theory. 



XXIIIrd International Symposium on Lattice Field Theory 

25-30 July 2005 

Trinity College, Dublin, Ireland 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http: //pos . sissa . it / 



Center Vortex Model for SU(3) Yang-Mills Theory 



Markus Quandt 



1. Introduction 

The vortex picture of the Yang-Mills vacuum was initially proposed as a possible mechanism 
of colour confinement. It is based on the idea that a random distribution of vortex colour flux is 
sufficient to effect an area law behaviour for large Wilson loops. Despite its simplicity and early 
successes [jl]], a precise confirmation of this scenario has been elusive for a long time. Recently, the 
picture has attracted renewed attention, mainly due to the advent of gauge fixing techniques which 
allow to detect center vortex structures directly within lattice configurations. Large-scale computer 
simulations have revealed ample evidence that the long-range properties of Yang-Mills theory can 
be accounted for in terms of vortices [0]. 

To complement the lattice approach, a random vortex world-surface model was introduced as 
an effective low-energy description of SU (2) Yang-Mills theory ; it has recently been extended 
to the gauge groups SU(3) 0, Bh. The fundamental assumption is that the long-range structure 
of Yang-Mills theory is dominated by thick, weakly interacting tubes of center flux which trace 
out closed surfaces in space-time. No analytical approach is known for such a quantum string 
ensemble. As a consequence, we realise our model on a space-time lattice with a fixed spacing 
a to represent the transverse thickness of vortices. Short-distance structures smaller than a (or, 
equivalently, momenta larger than n/a) cannot be resolved in this approach. 

I will discuss the physical foundations of the center vortex model and some technical aspects 
in the next section. The main part of the talk is section 3 which presents a selection of results 
obtained in the model. In Section 4, 1 conclude with a brief summary and some general comments. 



2. Phase diagram and choice of parameters 

Vortices are closed lines of colour flux in three space dimensions and, correspondingly, closed 
world-surfaces in space-time. Their dimensionality is unique in that they can have a topologically 
stable linking with Wilson loops. The colour flux placed on the vortex is quantized such that each 
linking with a Wilson loop W^) (i.e. each intersection point of the closed loop ^ with the vortex 
surface) contributes a center element of SU(N) to (WCrf)). To describe such intersections, the 
vortex surface must be defined on a lattice which is dual to the one where Wilson loops (and gauge 
connections) live. This emphasises the dual character of our model. 

Technically, we create random vortex surfaces on the dual lattice by assigning triality 1 to 
the elementary squares. The fixed lattice spacing a defines the minimal distance at which two 
intersection points can be resolved, i.e. the thickness of the vortex tube or sheet. Random surfaces 
created in this way are weighted by a model action containing a Nambu-Goto and a curvature term 
[Q], symbolically 




(2.1) 



The parameters are determined by measuring the zero-temperature string tension Oq (in units of 
the lattice spacing). One observes a deconfined region at large couplings (e,c) with a shallow 

'For the case G = SU(3), the center comprises three elements Z(3) = {1, e 2 ' 1 ''' 3 , e 4jr '/ 3 }. They are usually 
parametrised as z = exp(2jri/3 • q) in terms of the triality q £ {0, 1,2} defined modulo 3. The triality of an elemen- 
tary square determines the center element which a Wilson loop receives when intersecting the square. 
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Figure 1: String tension between two static colour charges (crosses) and spatial string tension (squares) 
as a function of temperature. Measurements were taken on a 16 3 x No lattice for the physical choice of 
parameters eq. (2.2). 



cross-over to a confined region at small (£,c) [|]|. In the latter domain, one can also effect a 
deconfinement phase transition by increasing the temperature of the simulation (i.e. by reducing the 
temporal lattice extension). Using an interpolation procedure, the ratio of the critical temperature, 
T c , to the zero-temperature string tension, Co, can be extracted at all couplings. Comparision with 
the known value T c /^/gq = 0.63 from St/ (3) lattice gauge theory then yields a line of possible 
physical values in parameter (e, c) space. As it turns out, the long-range physics are almost constant 
on this parameter line so that we are free to make the arbitrary choice, 

£ = 0, c = 0.21. (2.2) 

From the string tension Co a 2 measured at this point, the overall scale can be determined using the 
phenomenological value On = (440 MeV) 2 from QCD. This gives a vortex thickness of a = 0.39fm. 



3. Applications 

3.1 Spatial string tension 

The zero-temperature string-tension Co used to fix the parameters is extracted from Wilson 
loops extending in one space and the time direction. In addition, purely spatial Wilson loops 
can also be measured, although they do not have an immediate interpretation as an inter-quark 
potential. From lattice gauge theory, it is known that spatial Wilson loops exhibit an area law at 
all temperatures and the spatial string tension extracted from them persists in the high-temperature 
phase As shown in fig. [I], our model reproduces this strongly correlated hot phase; in fact, our 
results for the spatial string tension agree with ref. |§] to within 5%. This success can be attributed 
to the percolation behaviour of vortex clusters. 

3.2 Finite temperature phase transition and vortex branching 

The confinement mechanism is most clearly seen in the details of the phase transition at non- 
zero temperatures. Histograms of the action density measured at the critical temperature reveal a 
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Figure 2: The right-hand panel displays a four-dimensional lattice view, where a vortex surface, consisting 
of elementary squares, branches along av = 3 bond. The left panel shows the same situation projected to a 
3D slice. 

qualitative difference between the gauge groups SU (2) and SU (3): The SU (3) transition exhibits 
the shallow double-peak structure characteristic for a weak first order transition, while the SU (2) 
transition is continuous (second order) [|J]. This pattern is also seen in full lattice gauge theory. 

To understand the origin of this difference, we have to look at the vortex geometry in more 
detail. In D = 4 dimensions, each bond on the dual lattice is attached to six elementary squares, 
which are assigned trialities q 6 {0, 1,2} in our model. The local geometry of a world surface is 
then characterised by the number v = 0, ... ,6 of vortex surfaces (q ^ 0) meeting at a given bond. 
The odd values v = 3, 5 are not allowed in SU (2) and represent genuine SU (3) vortex branchings, 
cf . fig. ||. Note that the triality q = 1 , 2 of a vortex square can be reversed by flipping its orientation; 
moreover, triality is only conserved modulo 3 and a Dirac string (q = 3) may be added to any 
configuration. The branching in fig. ||| is thus equivalent to a situation where three q = 1 vortices 
emerge from a common point (or bond in D = 4). Branching lines in D = 4 are therefore similar to 
center monopole worldlines. 

The statistical distribution of branchings is best studied in 3D slices of the lattice, whence 
possible branching bonds are projected onto points of type v. From fig. ||, we conclude that the 
largest volume fraction in the confined phase corresponds to non-branching vortex matter (v = 2), 
with a considerable probabitlity of both self-intersection (v = 4, 6) and genuine branchings (v = 
3,5). Only 15% of the volume is not occupied by vortices (v = 0). In the deconfmed phase 
(T > T c ), the situation is qualitatively unchanged for time-slices, while space-slices show virtually 
zero branchings above T c . This can be understood if the vortices undergo a (de)percolation phase 
transition for T > T c and most vortex clusters wind directly around the short time direction [0]. 
The same phenomenon explains the persistence of the spatial string tension discussed above. 

3.3 Vortex free energy and the 't Hooft loop 

The 't Hooft loop can be viewed as a creation operator for (quantized) flux along a closed 
spatial line c to. It was formally introduced by 't Hooft in 1978 as a dual order parameter for the 
deconfinement transition on a space-time torus. Recently, explicit realisations of this construction 
have been given both in the continuum [|7|] and on the lattice [g], where maximally extended 't Hooft 
loops implement twisted boundary conditions. 

2 The case v = 1 describes a single vortex surface ending in the given bond; this is forbidden by Bianchi's identity. 
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Figure 3: Volume fractions occupied by points of a certain branching type v within 3D lattice slices. The 
left panel shows the distribution at zero temperature in the confined phase. The middle and right panel both 
correspond to T > T c , with the middle referring to a time slice and the right to a space slice. 

In the center vortex model, the 't Hooft loop is literally an (open) vortex creation operator 
[§] : Its action is to add a fixed triality q = 1 , 2 to each elementary square in a world-sheet over 
c £\ since triality is additive, this simply injects a center vortex of type q = 1,2 in the current 
configuration. The exact form of the world-sheet over ^€ is irrelevant (and center-gauge dependent), 
but for simplicity, we restrict ourselves to planar loops ^ and minimal surface sheets over them. 

The action penalty AS incurred by the vortex creation is related to the vortex free energy F 
via e~ F = (e~ AS ). As mentioned above, the free energy is expected to furnish an order parameter 
for the deconfinement transition with essentially the opposite behaviour as the Wilson loop. This 
is nicely confirmed in our model: The left panel of fig. ^| exhibits a linear rise of the free energy 
with the area of the 't Hooft loop at T > T c , 3 which allows to define a dual string tension Od in the 
deconfined phase. 

As we approach the phase transition from above, the dual string tension quickly vanishes 
(cf. right panel of fig. In the confined phase, the subleading perimeter law is hidden in the 
statistical noise, and the vortex free energy is consistent with zero (vortex condensation). Precise 
measurements close to the transition reveal a small discontinuity 

y/oB\ T=T ^ = (34.5 ±4.9) MeV (3. 1) 

in the free energy, which should be compared to the ordinary zero-temperature string tension On = 
(440MeV) 2 setting the overall scale. This demonstrates again the weakness of the first order 
transition for the colour group G = SU (3). Quantitatively, our findings are in fair agreement with 
lattice caluclations [Js[] which seem to favour a slightly larger value Od — (48 MeV) 2 . 

4. Further remarks and conclusions 

In this talk I have presented an effective model for the infra-red sector of SU (3) Yang Mills 
theory, based on random world surfaces carrying center flux. Assuming only that vortex sheets 
have a surface tension and stiffness, the model reproduces many non-trivial properties of long- 
range Yang-Mills theory. Among the examples discussed here are the spatial string tension, the 
order and strength of the phase transition for various gauge groups, and the discontinuity of the 
vortex free energy across the SU (3) phase transition. 

3 The systematic deviations can be understood in terms of Z(3) monopole correlations, see n for details 
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Figure 4: Left Panel: Free energy of 't Hooft loops in the deconfined phase (T/T c = 1.093), as a function 
of the minimal area over the loop. Right panel: The dual string tension Ou as a function of the temperature. 
Measurements were performed on a large 30 3 x iVo lattice with No = 1,2. 

The geometrical structure of vortex branchings is the key property in establishing first order 
behaviour for the SU (3) model. Geometry alone, however, is not sufficient to determine the physics 
of a vortex model; the effective vortex dynamics specific to each gauge group play a crucial role. 
In principle, the vortex dynamics would be determined by integrating out all non- vortex degrees of 
freedom; our project is essentially the reverse approach of guessing a model action and comparing 
to low-energy Yang-Mills theory. 



For SU(3), the results indicate that the simple two-operator action, eq. (2.1), is sufficient. 
With increasing complexity of the colour group, however, additional terms may arise which favour 
vortex branchings (center monopoles) explicitly; such terms tend to enhance the strength of the 
phase transition. In fact, preliminary investigations for G = SU (4) indicate that such terms are 
indeed necessary and may become even more pronounced with increasing number of colours. 
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